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We introduce for a general correlation scenario a new simulation model, a local quasi 
hidden variable (LqHV) model, where locality and the measure-theoretic structure 
inherent to an LHV model are preserved but positivity of a simulation measure 
is dropped. We specify a necessary and sufficient condition for LqHV modelling 
and, based on this, prove that every quantum correlation scenario admits an LqHV 
simulation. Via the LqHV approach, we construct analogs of Bell-type inequalities for 
an iV-partite quantum state and find a new analytical upper bound on the maximal 
violation by an iV-partite quantum state of Si x • • • x Sjv-setting Bell-type inequalities 
- either on correlation functions or on joint probabilities and for outcomes of an 
arbitrary spectral type, discrete or continuous. This general analytical upper bound 
is expressed in terms of the new state dilation characteristics introduced in the present 
paper and not only traces quantum states admitting an Si x ■ ■ • x SV-setting LHV 
description but also leads to the new exact numerical upper estimates on the maximal 
Bell violations for concrete iV-partite quantum states used in quantum information 
processing and for an arbitrary iV-partite quantum state. We, in particular, prove 
that violation by an iV-partite quantum state of an arbitrary Bell-type inequality 
(either on correlation functions or on joint probabilities) for S settings per site cannot 
exceed (25'— l)^ -1 even in case of an infinite dimensional quantum state and infinitely 
many outcomes. 
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I. INTRODUCTION 

The seminal papers of Einstein, Podolsky and Rosen 1 (EPR) and Bel^ are still ones 
of most cited in quantum information. In Ref.-, Einstein, Podolsky and Rosen argued that 
locality of measurements performed by spatially separated parties on perfectly correlated 
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quantum events implies the "simultaneous reality - and thus definite values" of physical 
quantities described by noncommuting quantum observables. Based on this argument con- 
tradicting, however, the quantum formalism and referred to as the EPR paradox, Einstein, 
Podolsky and Rosen expressed in Ref.- their belief on a possibility of a hidden variable 
account of quantum measurements. 

Analyzing this EPR belief in 1964 - 1966, Bell explicitly constructed^ the hidden variable 
(HV) model reproducing the statistical properties of all quantum observables for a qubit 
but, however, proved 3 that, for bipartite measurements on a two-qubit system in the singlet 
state, a local hidden variable description (LHV) disagrees with the statistical predictions 
of quantum theory. Based on these results, Bell concluded^ that the EPR paradox should 
be resolved specifically via violation of locality under bipartite quantum measurements and 
that " ...non- locality is deeply rooted in quantum mechanics itself and will persist in any 
completion" . 

Ever since 1964, the conceptual and mathematical aspects of the probabilistic description 
of multipartite quantum measurements have been analyzed in a plenty of papers, see, for 
example, articles^— and references therein. Nevertheless, as it has been recently noted by 
Gisin 11 , in this field there are still "many questions, a few answers". 

It was, for example, proved by Werner 8 that there exist finite dimensional nonsepara- 
ble bipartite quantum states admitting an LHV description under all projective bipartite 
quantum measurements with an arbitrary number of measurement settings at each site. It 
was also shown in Refs.— ^2r— that some nonseparable bipartite quantum states admit an 
LHV description only under correlation scenarios with specific numbers of measurements 
at iV sites. However, until now it is not still known what state parameter quantitatively 
determines violation by an iV-partite quantum state of Bell-type inequalities 16 - constraints 
specifying scenarios admitting an LHV description and named after the seminal result^ of 
Bell. 

Nowadays, it is also clear— that though multipartite quantum measurements do not need 
to be local in the sense of Bell, they are, however, local in the sense meant by Einstein et al 
in Ref™. The difference between the general nonsignaling condition, the EPR locality and 
Bell's locality is analyzed in Ref. 10 . Thus, the term "a nonlocal quantum state" widespread 
in quantum information means now only that this state does not admit an LHV description 
and, therefore, violates some Bell-type inequality. 
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This takes us back to the EPR locality argument^ and asks - if it is possible to construct 
for a quantum correlation scenario a simulation model which would be (i) local in the sense 
meant by Einstein, Podolsky and Rosen^; (ii) similar by its measure-theoretic construction to 
the concept of an LHV model and (iii) incorporate the latter only as a particular case. This 
problem is also urgent for all multipartite correlation scenarios (not necessarily quantum) 
specified not in terms of a single probability space. The latter is one of the main notions of 
the conventional probability theory. 

Apart from the purely theoretical interest, such a local simulation model could also single 
out a state parameter characterizing quantitatively violations of Bell-type inequalities by a 
multipartite quantum state - the problem discussed in the literature ever since the seminal 
result of Tsirelson^. 

Note that though, for correlation bipartite Bell-type inequalities, quantum violations are 
upper bounded- by the Grothendieck's constant^ independently on a dimension of a bipar- 
tite quantum state and numbers of settings and outcomes at each site, this is not already the 
case for bipartite Bell-type inequalities on joint probabilities. Since Bell-type inequalities 
are now widely used in many quantum information tasks^ 1 ^, bounds on quantum viola- 
tions of Bell-type inequalities have been recently intensively discussed in the literature both 
computationally— and theoretically^— and it has been founds that some tripartite quan- 
tum states "can lead to arbitrarily large" violations of correlation Bell-type inequalities. For 
an iV-partite quantum state, bounds on violation of a Bell-type inequality of an arbitrary 
type (either on correlation functions or on joint probabilities) have not been reported^ in 
the literature. 

In the present paper, we introduce for the probabilistic description of a general correlation 
scenario a new simulation model, a local quasi hidden variable (LqHV) model, where locality 
and the measure-theoretic structure inherent to an LHV model are preserved but positivity 
of a simulation measure is dropped. We prove that every quantum correlation scenario 
admits the simulation in LqHV terms and construct via the LqHV approach analogs of Bell- 
type inequalities for an iV-partite quantum state. This allows us to find the new analytical 
and numerical upper bounds on the maximal violation by an iV-partite quantum state of 
all Bell-type inequalities - either on correlation functions or on joint probabilities and for 
outcomes of an arbitrary spectral type, discrete or continuous. The paper is organized as 
follows. 
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In section 2, for our consideration in sections 5, 6, we specify some new dilation charac- 
teristics of an iV-partite quantum state and discuss their properties. 

In section 3, we introduce for a general iV-partite correlation scenario with S n measure- 
ments at each n-th site the notion of an LqHV model and specify a necessary and sufficient 
condition for LqHV modelling. 

In section 4, we recall for an Si x • • • x Sjv-setting correlation scenario with outcomes 
of an arbitrary type, discrete or continuous, the general form of all Bell-type inequalities - 
either on correlation functions or on joint probabilities. 

In section 5, we prove that every quantum Si x • • • x S^r-setting correlation scenario 
admits an LqHV model and introduce, for an iV-partite quantum state, the exact analytical 
upper bound on the state parameter specifying a possibility of its Si x • • • x Sjv-setting LHV 
description. 

In section 6, via the LqHV approach, we construct analogs of Bell-type inequalities for 
an iV-partite quantum state and find the new analytical and numerical upper bounds on 
the maximal violation by an iV-partite quantum state of all S\ X • • • x S^r-setting Bell- 
type inequalities. The comparison of our exact general iV-partite numerical upper estimate 
specified for N = 2, 3 with the bipartite and tripartite numerical estimates reported in the 
literature is given in section 6.2. 

In section 7, we summarize the main results of the present paper. 

In appendices A, B, C, we present proofs of some statements formulated in sections 2, 5 
and 6, respectively. 

II. PRELIMINARIES: SOURCE OPERATORS, TENSOR POSITIVITY, 
THE COVERING NORM 

In this section, for our consideration in sections 5, 6, we specify the notion of a source 
operator—^ for an iV-partite state, the notion of tensor positivity^ and introduce a new 
norm, the covering norm, on the space of all self-adjoint trace class operators on a tensor 
product Hilbert space. 

For a quantum state p on a complex separable Hilbert space H\ ® • • • ® and arbitrary 
positive integers S\,...,Sn > 1, denote by T^... xg a self-adjoint trace class operator, 
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defined on T-L\ 



T-if i SN and satisfying the relation 



tr 



j'(p) 



= tr 

fc 1 = 0,...,(S' 1 -l),...,A;j V = 



(1) 



for all bounded linear operators Xi,...,Xn on Hilbert spaces "Hi, Hn, respectively. In 



I]), we set I n m ®X n \ k=0 = X n ®l u m \ k=0 := X n . Clearly, tr[T^ 



rtn 



' Si x-xSjy J 



1 and TOO 



Definition 1 (Source operators For a state p on a Hilbert spaced 7i 



and arbitrary positive integers Si, > 1, we call each of self- adjoint trace class operators 



TW „ on ftf 



<E> ■ • ■ <S> T-L% Sn satisfying relation (QP as an Si x ■ • • x S^-setting source 



operator for state p. 

For a source operator T, its the trace norm 

UTI^ := tr[|T|] = tr [T + +T~]=1 + 2tr[T~] > 1. (2) 

Here, T± > 0, T+T~ = T~T+ = are positive trace class operators in the spectral decom- 
position T = T + — T~ and |T| := a/T^ = T + + T~. 



Proposition 1 For even/ state p on a Hilbert space Hi 
integers Si,...Sn > 1, t/^ere exists a source operator T,W 



"Htv fln< ^ arbitrary positive 



Proof. For a bipartite case, this statement has been proved for settings 1 x 2 , 2 x 1 by 
proposition 1 in Ref.—. This proof was further generalized in appendix of Ref.— for arbitrary 
1 x S2, Si x 1. The proof for a general iV-partite case with setting Si x ... x Sn is presented 
in appendix A. ■ 

If „ is a source operator for state p, then each of its reduced ( „ ) on a 

SiX-xSjv V SlX "' xS W red 

Hilbert space Uf Ll g> • • - 
source operator for state p and 



T-L% nN , with 1 < L n < S" n , constitutes an Li x ■•■ x Ljv-setting 



1 < 



( T (p) 



< 



j"(p) 

Si x-'-xSjv 



(3) 



In order to analyze situations where, for a source operator T, relation tr[T{Xi <g> ■ ■ ■ 
<g)X m }] > holds for arbitrary positive operators X\, ...,X m , we specify the following general 
notion. 
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Definition 2 (Tensor positivity—) We call a bounded linear operator Z on a Hilbert 
space Qi ® • ■ • ® Q m , m> 1, as tensor positive and denote it by Z > if the scalar product 

® ■ ■ ■ ® i> m , z i>i ® • • • ® > o (4) 

/or arbitrary ip x E Qi, ...,ip m G £ m . 

Remark 1 For space Q\ ® Q2, the notion of tensor positivity is similar by its meaning to 
"block-positivity" in Ref.—. We, however, consider that, for a tensor product of any number 
of arbitrary Hilbert spaces, possibly infinite dimensional, our term "tensor positivity" is more 
suitable. 

For m = l, tensor positivity is equivalent to positivity. For m > 2, positivity implies 
tensor positivity but not vice versa. For example, operator V(ip 1 <g> ip 2 ) := ^2 ® ^1 011 space 
W <8> W is tensor positive but not positive. 

Since on a complex separable Hilbert space, every positive operator is self-adjoint, from 
(J3J) and the spectral theorem it follows that, for a trace class tensor positive operator W > 
on a complex separable Hilbert space Qi®- ■ -®Gmi relation tr [W{Xl ® ■ ■ ■ ® X m }] > holds 

for arbitrary positive operators X 1 ,...,X m on spaces Qi,..,Q m , respectively. In particular, 

® 

tr[VK] > 0, for each trace class W > 0. 

If a trace class operator on Q\ ® ■ • • ® Q m is tensor positive, then any of its reduced 
operators is also tensor positive. The converse of this statement is not true. 

Coming back to source operators, we stress that though, for every iV-partite state, a 
source operator exists (see proposition 1) for every setting S\ x • • • x Sn, an arbitrary N- 
partite state does not need to have a tensor positive source operator^. 

For example, every separable iV-partite state p = Yl a iPi ® • ■ ■ ® P% cti > 0, Yl a>i = 1, 
has a positive source operator 

E^fV" 8 ^) 8 *" ( 5 ) 

for arbitrary Si, > 1. However, a nonseparable state does not need to have a tensor 
positive source operator even for at least one setting. In Refs.— we present examples 
of source operators for some nonseparable bipartite states and single out the state parameters 
for which these source operators become tensor positive. 

Suppose now that we want to decompose a source operator into two tensor positive 
operators. The following new notion allows us to consider such decompositions. 
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Definition 3 (Coverings) For a self- adjoint bounded linear operator Z on a Hilbert space 
Gi ® • ■ ■ ® Qm, Tn > 1, we call a tensor positive operator Z cov on Q\ <S> ■ ■ • <8> Q m satisfying 
relations 

Z cov ±Z>0 (6) 

as a covering of Z. 

If Z is tensor positive, then it, itself, represents one of its coverings. In view of fl6]), every 
self-adjoint bounded linear operator Z on Q\ £g> ■ • • <g> Q m admits the decomposition 

Z= l - (Z cov + Z)- 1 - (Z cov - Z) (7) 

via tensor positive operators Z cov ± Z > 0, where Z cov is any of its coverings. 

For a source operator, we are interested in its trace class coverings. Denote by 7g 1( g>...®g m 
the linear space of all trace class operators on a Hilbert space Qi®- ■ -®^ m and by Tg^... l ^g m C 
T<3i®i"®Qm ~ the subspace of all self-adjoint trace class operators. 

Proposition 2 For every self-adjoint trace class operator W on a Hilbert space Q\®- ■ -®Q m , 
there exists a trace class covering W cov . 

Proof. For an operator W G Tg%...^,g m i consider its spectral decomposition W = W + — W~ , 
where W ± > 0, W + W~ = W~W + = 0. The positive trace class operator \W\ := VW^ = 
W + + W~, with tr[|W|] := < oo, constitutes a trace class covering of W. This proves 

the statement. ■ 

From it follows that, for every self-adjoint trace class operator W, the relation 

trfWcJ > \ti[W]\ >0 (8) 

holds for each of its trace class coverings W cov . 

Thus, for every W G 7g'g!, l8 g , the set {VFco,; G Tg l( ^...^g m } 3 \W\ contains at least one 
element and tr [Wco„] > for each covering W cov . Therefore, we can introduce on space 
Tgi®-®g m the following function 

f ^ W) '-=w ,i nf *riW cov ]>0, VWgTS..^> (9) 



s 



and relations 



f(W)=0 W = 0, (10) 

fiptW) = \a\ f(W), Va G R, 
f(W 1 + W 2 )<f(W 1 ) + f(W 2 ), 

hold for all W, Wi, W 2 G 7g^...^g m - The first of these relations follows from property 1 in 
lemma 1 below. The second relation - from (0), fl9]). For the proof of the last relation 
in (IT01) . we note that, for arbitrary trace class coverings (Wi) cov , (W 2 )cov of operators Wi, 
^ 2 e ^Qt®-®Q ' °P era t° r (Wi)ow + (W^co-u constitutes a possible trace class covering of 
VFi + W 2 . Hence, {(Wi + W 2 ) cov } D {(Wi) cov + (W 2 ) cov } and, taking this inclusion into the 
account in infimum specifying f{W\ + W 2 ), we come to the third relation in (HTJ]) . 
In view of relations (ITU]) , function OH]) constitutes a norm on space Tgl®...®g m 

Definition 4 VKe re/er to norm (G?P as i/ie covering norm and denote it by 



\W\ 



inf 



tr[W 



The following general properties of the covering norm are proved in Appendix A. 



(11) 



Lemma 1 For a self-adjoint trace class operator W on a Hilbert space Q\ ( 

(1) MW]\ <sup|tr[W{Xi® ...®X m }]| < ||W|| C(TO < HWlli, 

where supremum is taken over all self-adjoint bounded linear operators Xi, 
erator norms \\Xj\\ = 1 on spaces Q\, ...,Q m , respectively; 

(2) ifW>0, then \\W\\ cov = tr[W]; 



(3) \%x[W]\ < \\W red \ 



< 



for each operator W re d reduced from W. 



Qm ■ 



, X m with op- 



As an example, consider the self-adjoint operator V(ip 1 <g> ip 2 ) = ip 2 £g> ip 1 on C d <S> C d . For 
this operator, the covering norm ||V|| cow = d while the trace norm {{V^ = d 2 . 



For an S*! x ... x SV-setting source operator I^J xS ^ for a state p on a Hilbert space 



•Hi 



"Hat, lemma 1 implies 



1 < 



t (p) > o 



j"0) 

Si x ■■■xS N 



J>{p) 

SiX'-'XSjv 



< 



Si x-'-xSjv 



(12) 
(13) 
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and 



1 < 



y(p) 

Si x---xS N 



< 

cov 



y(p) 



(14) 



for each source operator [Tff) ) reduced from a source operator 



' red 



III. LQHV MODELLING OF A GENERAL CORRELATION SCENARIO 

Consider an AT-partite correlation scenario, where each n-th of iV > 2 parties (players) 
performs S n > 1 measurements with outcomes A n G A n of an arbitrary type and F\ n is 
a cr-algebra of events _F ra C A n observed at n-th site. For the general framework on the 
probabilistic description of multipartite correlation scenarios, see Ref.— 

We label each measurement at n-th site by a positive integer s n = l,...,S n and each 
of TV-partite joint measurements, induced by this correlation scenario and with outcomes 
(Ai, . . . , Ajv) G Ai x • • • x Ajv - by an A^-tuple (si, sjv), where n-th component refers to 
a measurement at n-th site. 

For concreteness, we further refer to an Si x • • • x SV-setting correlation scenario with 
outcomes in Ai x • • • x An by symbol 

£ s ,a, S := Si x ■ ■ • x S N , A := Ax x ■ ■ ■ x Ajv, (15) 

if \ 

and denote by -P( Sl S A SJV ) a probability measure, describing an A^-partite joint measurement 
(si, sn) of scenario and defined on the direct product (Ai x---xAjv, J~a x ®- • ■^)J i 'a n ) of 
measurable spaces (A n , J-a„), n = 1, N. Recall^ that the product cr-algebra Ta^- • -^^An 
is the smallest cr-algebra generated by the set of all rectangles Fi x • • • x Fn C Aj X • ■ • x 
with measurable "sides" F n G .Fa,,, n = 1, N. 

In what follows, we consider only standard measurable spaces. In this case, each (A n , Fa„) 
is Borel isomorphic to a measurable space (X n ,Bx n ), where X n E Br is a Borel subset of R 
and Bx n '■= Br fl X n is the trace on X n of the Borel cr-algebra Br on R. 

For a general correlation scenario £s,a, let us introduce the following new type of simu- 
lation models. 



Definition 5 We say that an Si x ... x Sn -setting correlation scenario £s,a, with joint prob- 
ability distributions P^' A SN y si = 1, Si,...,sn = 1, Sn, and outcomes (Ai, . . . , Ajv) G 
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Ai x • • • x Atv := A, admits a local quasi hidden variable (LqHV) model if all of its joint 
probability distributions admit the representation 



x-..xF N ) = J P^iF^u) ■ . . . ■ P^\f n \uj) v Ss Jdu), 



(16) 



n 

^1 £ ^Ai, •••> F N e F^ 



in terms of a single measure space [Q, Fn, ^£ sa ) , itf^/i a normalized bounded real-valued 
measure Vs Sh , and conditional probability measures Pn n \-\uj) : F\ n — > [0, 1], defined v Sgt ^- 
a.e. (almost everywhere) on Q and such that, for each s n = 1, S n and every n — 1, N, 
function Pn (F n \-) : O — > [0, 1] measurable for all F n G J r A n . 

Notation 1 In a triple (Q, Fn, v) representing a measure space, Q is a non-empty set, Fn is 
a o -algebra of subsets offl and v is a measure on a measurable space (Q, Fn) . A real-valued 
measure v is called normalized if u(Q) = 1 and bounded if \v{F)\ < M < oo for all F e .Fn- 



We stress that, in an LqHV model ([16]) . measure V£ SA has a simulation character and 
may, in general, depend (via the lower index £,s,a) on measurement settings at all (or some) 
sites, as an example, see measure f|40|) . 

From f)16p it follows that a correlation scenario £$,a admitting an LqHV model satisfies 
the general nonsignaling condition specified by definition 1 (Eq. (10)) in Ref.—. 

If, for a correlation scenario £s,a, there exists representation f |T6|) . where a normalized 
real- valued measure vz s K is positive and, hence, is a probability measure, then this scenario 
admits an LHV model formulated for a general case by definition 4 (Eq. (26)) in Ref.—. 

Remark 2 Recall that a bounded real-valued measure v on a measurable space (Q, Fn) 
admits the Jordan decomposition v = v + — v = via positive measures 

u + (F):= sup u(F'), u~(F):=- mi u(F), VF G Fn, (17) 

F'e^F'CF F'GFntF'CF 

with disjoint supports. The sum (z/ + (fi) + v~(Q)) coincides with the total variation \u\ (Q) 
of measure v on Q, which is defined by relation 

m 

sup^K^)|:=M(^HL,, (18) 

i=l 

where supremum is taken over all finite systems {Fj} of disjoint sets in Fn. For a bounded 
measure v, its total variation \\v\\ var < oo and \\-\\ var constitutes a norm, the total variation 
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norm, on the linear space of all bounded real-valued measures on a measurable space (O, J-"n). 
Thus, for a bounded real-valued measure z/, we have 



v+(n) + v-(n) = \\v\\ var . (19) 

// a bounded real-valued measure v is normalized, then 

IHI w = i + 2f-(n)>i. (20) 

A normalized bounded real-valued measure v is a probability measure iff \v\ var = 1- Note 
that relation 

sup HF)\<\\v\\ var <2 sup \u(F)\ (21) 
holds for every real-valued measure v. 

From the Jordan decomposition for measure ve SA it follows that if a correlation see- 



nario admits an LqHV model (|T6|) . then each of its joint probability distributions P^' K SN ) 
can be expressed via the affine combination of some LHV distributions Pr s ^' A SN \ that are 

>0n)/ 



represented in (jT6l) by the same conditional measures Pn n (-\uj). On the other hand, if a 
correlation scenario with a finite number of outcomes at each site admits the affine model in 
the sense of Ref.— , then this scenario admits the special LqHV model, where measure ug 3A 
is given by the affine combination of discrete probability measures and each P!f n) (F n \u) 
has the particular form x _ ± (u), F n G J-A n , where f niSn : Q — > A n is some measurable 
function, fnl n (Fn) := {w £ fn,s n (u) G F n } and Xd(') ^ s the indicator function of a subset 
DCfi, that is: Xd{ u ) = 1 if w G -D and Xd{ u ) = if u ^ D. 

Thus, an LqHV model incorporates as particular cases and generalizes in one whole both 
types of simulation models discussed in the literature - an LHV model and an affine model. 
Note that the latter model is, in principle, built up on the concept of an LHV model. 

We stress that, in an LqHV model, locality and the measure-theoretic structure inherent 
to an LHV model are preserved. 

The following general theorem introduces a necessary and sufficient condition for LqHV 
modelling. 



Theorem 1 An S\ x ... x Siy-setting correlation scenario £s,a admits an LqHV model (B 
if and only if, on the direct product space (Af 1 x • • ■ x 

^Ni^Ai 1 ® • • • ® ^a^), there exists 
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a normalized bounded real-valued measure^ 

fi £sA (d\? x • • • x d\[ Sl) x • • • x dA$ x • • • x dA$ w) ) , (22) 
X^ n) eA n , s ft = l, ...,£„, n = l,...,JV, 

returning all joint probability distributions P^' A Sjv ) 0/ scenario E$,a as the corresponding 
marginals. 

Proof. Let scenario £s,a admit an LqHV model ffTBI) Then the normalized real-valued 
measure 

I { J] P^\d\^\uj)}u £s ^) (23) 

** ~ 1 C 1 AT 



s n =l,...,S n , n=l,...,N 



on (Af 1 x • • • x A^pff 1 <g> • • • ® P? 5Ar ) returns all distributions P ( f s ' A) s , of scenario £ s a 
as the corresponding marginals. The total variation of measure (1231) is upper bounded by 
It^saII < 00 ' so that, in view of relation (12T]) . this measure is bounded. 

In order to prove the sufficiency part of theorem 1, let there exist a normalized bounded 
real- valued measure u Fo returning all probability distributions P, S,A of scenario as 
the corresponding marginals. This means that the representation 

P^L) (FiX---xF N ) = J Xfi (X^) ..... Xfn (X^) H £sA (d\? (24) 

x • • ■ x dX[ Sl) x • ■ ■ x dA$ x • • ■ x dA^), 

S\ = 1, Si, Sat = 1, Sat, 

holds for all Pi G Pa i; ...,P/v G Pa^- Representation ( 124)) constitutes a particular case of 
the LqHV representation ([TBI) specified with 

u/ = ^A[ \ x[ x \ A^, A Ar JV ' ) J , (25) 
= Af 1 x ... x A^, P n , = Pff 1 ® • • • ® Pf^, 

This proves the statement. ■ 

The following corollary of theorem 1 corresponds to the statements (a), (c) of the general 
theorem 1 on LHV modelling in Ref.— . 
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Corollary 1 An Si X ... X Sn -setting correlation scenario £$,a admits an LHV model if and 
only if there exists a probability measure a4 s A on space (Af 1 x • • ■ x A^, -Ff-f 1 ® • • • ® Ff N N ) 
returning all joint probability distributions P^' A Sjv ) of scenario £$,a a s the corresponding 
marginals. 

Proof. If scenario £s,a admits an LHV model, then there exists representation fj!6j) with 
some probability measure is' and, for this is', the constructed normalized measure (I23p is 
a probability one. Conversely, let there exist a probability measure /i^ SA returning all 

( F \ 

distributions P^' A Sjv ) of scenario £$,a as the corresponding marginals. Then representation 
( 124)) with probability measure /jf £ A constitutes a particular LHV model. ■ 

IV. BELL-TYPE INEQUALITIES 

For an Si x • • • x ^-setting correlation scenario £s,a, with joint probability distributions 
P(si sjv) anc ^ outcomes (Ai, ■••) Ajv) G Aj x • • • x Atv := A, consider a linear combination 

E (^,..^)(Ai,..-,Ai,)) (26) 

s u ...,s N S ' A 

of averages 

(V'(. 1> ...,. jr }(A 1 ,.-,A JV )) £flA (27) 
:= / ^( Sll ..., SiV )(Ai, Av)P£i Jv) (dA 1 x ... x dAjv) 

A 

arising under joint measurements (si, . . . , s N ) and specified by a family 

*5,A := V0( Sl ,..., Sjv ), si = 1,..., 51,..., sjy = 1,...,Sjv} (28) 

of bounded measurable real- valued functions ip( slt ^ iS )(Ai, Ajv) on the direct product (A x x 
• - - x Ajv, J 7 a 1 <8> • • ■ <8> of measurable spaces (A n , Ja„)> n — 1, V. 

If, in (127)1 . function ^(si,...^ ) nas t ne product form sl (Ai) - ... -0 S (Ajv), then, depending 
on a concrete choice of functions 4> Sn (X n ), for a joint measurement (si, s N ), the average 

/ sl (Ai) ■ - ■ 4> Sn (X n ) P^ SN) (d\i x ... x dAjv) (29) 

A 
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may refer to either the joint probability 



(30) 



of events Fi G J-a 17 Fn G J-\ n observed at the corresponding sites or if outcomes are 
real-valued and bounded - to the expectation value 



( s »i) \( S71 m) 

At 



£s,a 



I An 1 -...-A nM P£^ ) (dA 1 x.--xdA iV ) 



(31) 



of the product A ni - ... - A n of outcomes observed at arbitrary M < N sites 1 < n\ < ... < 
n-M < N. For M > 2, the expectation value (1311) is referred to (in quantum information) as 
a correlation function. A correlation function for an iV-partite joint measurement is called 
full if. in (I5TI). M= N. 

If a correlation scenario £$ t A admits an LHV model, then every linear combination (126]) 
of averages satisfies the tight LHV constraints 17 



< A < E (^,.,s n) (^..,Xn)) p ^ \ lhv <B 



Sl,...,S. 



7SUp 



(32) 



where the LHV constants £>^ A and <BJ^ A constitute, correspondingly, supremum and infi- 
mum of fl26|) over all LHV scenarios and have the form: 



^*S P A :=SU P E \^(^,...^ iv )( A l'--' A ^' 



flhv 

°S,A Si 



flhv 



sup E VW.^)^"-'^)' 



Al Sn) 6A„,Vs„,Vn Sl ,...,s 



(33) 



15, 



inf 



M E (^ 1 ,...,s Ar )( A i,-,v 



flhv 
C S,A 



A™eA„,Vs„,Vn 



51, ...,s» 



Constraints (1321) imply 



E (^i,.., Sjv )( A i'-' v 



31,. 



LHV 



(34) 



where 



15 



By SA ■— max 



(35) 



= sup 

Ai Sn) GA n ,Vs„,Vn 

Inequalities f l32]) have been introduced in Ref.— and represent the general form of all 
unconditional 37 tight linear LHV constraints on correlation functions and joint probabilities 
for an 5*i x • • ■ x SV-setting correlation scenario with outcomes of an arbitrary type, discrete 
or continuous. 

Note that some of the LHV constraints (I3"21) may be fulfilled for a wider (than LHV) class 
of correlation scenarios. This is, for example, the case for those LHV constraints on joint 
probabilities that follow explicitly from positivity and nonsignaling of probability distribu- 
tions Pfc*; aif) and are, therefore, fulfilled for any nonsignaling scenario £s,a- Moreover, for 
some *&s,A, the corresponding constraints ( |32|) may be simply trivial - in the sense that these 
constraints are fulfilled for each scenario £$,A- For example, if we specify fl32|) with functions 
■0(si,...,s JV )( / ^i> •••> ^n) = 1) V(Ai, Ajv) G A, for all joint measurements (si, . . . , s N ), then 

41= £ x ")) es r^l =Si --- SN (36) 

si,...,s N 

holds for every scenario 

If, however, an LHV constraint may be violated in a non-LHV case, then it is generally 
named after Bell due to his seminal result in Ref.-. 

Definition 6 Each of the tight linear LHV constraints HM) that may be violated under a 
non-LHV correlation scenario is referred to as a Bell-type (equivalently, Bell) inequality. 

As it is discussed in section 3 of Ref.—, the general form (1321) covers in a unified manner 
all unconditional Bell-type inequalities that were introduced via a variety of methods ever 
since the seminal publication of Bell 3 ". Note that the original Bell inequality 3 *, discussed 
recently in Ref.—, constitutes an example of conditional Bell-type inequalities. 

V. LQHV MODELLING OF A QUANTUM CORRELATION SCENARIO 

Let, under an S\ X • • • x S^r-setting correlation scenario, each A^-partite joint measurement 
(si, sn) be performed on a quantum state p on a Hilbert space 1-Li®- ■ -^Hn and described 

16 



s u ...,s N 



by joint probability measures 

tr[p{MS Sl) (dA!) ® • • • <g> M^dA^)}] (37) 

on the measurable space (Aj x ■ ■ • x Aat, ^Ai ® • • ■ ^^a^)- Here, each M^"' 1 is a normalized 
positive operator-valued (POV) measure on a measurable space (A n , Fa^) representing on 

(s ) 

a Hilbert space 7i n a quantum measurement s n at n-th site. For a POV measure 

MF, all 

its values M„ (P n ), F n G Ja„, are positive operators on H„ and (A n ) = I^ n . 

We specify this quantum Si x • • • x S^-setting correlation scenario by symbol £ p ,m sa , 
where 

M 5 ,a := {M^\ s n = 1, .., S n , n = 1, iV} (38) 
is a collection of POV measures describing this quantum scenario and denote by 

P £ M ,^)VAi x ■ ■ ■ x dA*) := tr [p {M^idX,) ® ■ ■ • ® M^dA*)}] , (39) 

Si = 1, Si, Sat = 1, SjV, 

its joint probability distributions (1571) . 

Theorem 2 For even/ state p on a Hilbert space Hi ® • • • ® "Hjv arbitrary positive 
integers Si, Sat > 1 ; eac/i quantum Si x ... x SN-setting correlation scenario £ p ,m sa , with 
joint probability distributions (E5j) and outcomes of an arbitrary spectral type, discrete or 
continuous, admits an LqHV model. 

Proof. For a state p on Hi <8> ■ ■ ■ <S> Hat, let Tg^ x _ xS be an Si x ... x SAr-setting source 
operator on space H,f Sl <8> • • ■®'H% Sn , see definition 1 in section 2. For each scenario ^ Pi m sa , 
the normalized real-valued measure 

^,M S , A ) / dA (l) x . ^ x dA (5 x _ x ^(1) x . . . x dX (S N )\ (4Q) 
T S 1 x---xS N V / 

== tr[T|f x ... x5iV {MS 1} (dAS 1} ) ® • • • ® N^dA? ) 
® • • - ® M$(dA$) ® • • • ® M^CdA^ 5 )}] 

on the direct product space (Af 1 x ■ • • x A^, Jff 1 <g> ■ ■ ■ ® J 7 ^) returns all joint probability 
distributions P, ^\ of scenario £ P! m s ,a as the corresponding marginals. Due to bound 
(IBlOp proved in appendix B and relation (fT2"j) . the total variation norm of measure PHI) is 

17 



upper bounded by | \T^ x ... xSn \\i < oo. This and relation (|2ip imply that the normalized real- 
valued measure ( 1301) is bounded. Thus, for each quantum scenario £ Pi m sa > the constructed 
measure p^ P i^ s ' A ^ satisfies the sufficiency condition of theorem 1 on LqHV modelling. This 
proves the statement. ■ 

If, for a state p, every quantum scenario £ p ,m sa (i-e. for an arbitrary collection M^a of 
POV measures and an arbitrary outcome set A) admits an LHV model, then, according to 
our terminology in Ref.— , this state p admits the Si x • ■ • x S^-setting LHV description. 
In the latter case, state p admits^ an Li x - • • x Ljy-setting LHV description for all L\ < 
Si,...,Ln < Sn, but does not need to admit the LHV description whenever at least one 
L n > S n . 

Via a similar terminology for the LqHV case, theorem 2 reads - every N -partite quantum 
state p admits an S\ x • • • x Siy-setting LqHV description for arbitrary numbers Si, Sn 
of measurements at N sites. 

In view of theorems 1, 2, corollary 1 and relation [201 let us introduce, for a quantum 
correlation scenario £ p ,m sa , the parameter 



7<? „ := inf 



a. A 



> 1, (41) 



var 



where, infimum is taken over all normalized bounded real-valued measures u r , each 
returning all distributions P^' of scenario £ p ,m s a as the corresponding marginals. 
The following lemma is proved in appendix B. 

Lemma 2 A quantum correlation scenario £ p ,m S a admits an LHV model if and only if 

7(7 = 1. 

Introduce also the state parameters 

^ ) .. x5jv :=sup7, p , MsA >l, (42) 
^ x ... X 5 w :=supT^.. xSw >l. (43) 

A 

Proposition 3 (a) For a state p on Hi <g> • • • (g) Hn, each quantum Li x • • • x L^-setting 
scenario £ p ,m l a , with L x < Si, L N < Sn and an outcome set A, admits an LHV model if 
and only ifT^.. x3N = 1; 

(b) A state p on Hi <8> • • • ®"Hjv admits an Si x • • • x Sn -setting LHV description if and only 
if SiX-xS N = 1- 
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Proof. If each scenario £ p ,m l a , where Li < Si, L N < Sn, admits an LHV model, then, by 
lemma 2, 7^ M ^ = 1 for all collections M Lj a of POV measures, where L x < S±, ...,L N < Sn- 
Hence, due to its definition fj42]) . parameter T<£' x ,,, xS = 1. Conversely, let T<f' x ,,, xS = 1. 
Then, in view of fj42j) . ( 14T]) . 7^ M =1 for all collections M^a of POV measures. By lemma 
2, this implies LHV modelling of every quantum correlation scenario £ Pi m S a- By proposition 
3 in Ref.— , the latter, in turn, implies LHV modelling of each quantum correlation scenario 
£p,M LA whh settings L\ < Si,...,L N < Sn- This proves the sufficiency part of statement 

(a) . Statement (b) is proved quite similarly. " ■ 

Thus, for an iV-partite state p, it is specifically the state parameter T^ x ... x5 that 
determines quantitatively a possibility of an LHV description of all quantum scenarios fl39|) 
with settings up to setting 5iX---x Sn and outcomes of an arbitrary type. 

Proposition 4 For a state p on a Hilbert space Hi® - ■ -®1-Ln and arbitrary positive integers 
Si, Sn > 1, 

, < T (p,A) < T (p) < . r nrpip) II (aa\ 

1 — L S 1 x-xS N — L S!X-xS N — - . 1111 \\ 1 S 1 x-xlx-xS N ll«w" K^J 

T S 1 x---xlx---xS N ' ^ n X 
t 

n 

where (i) infimum is taken over all source operators Tg^x— xix— xsw w ^ on ^V one setting 

t 

at n-t/i szie and over all n = 1,...,N; (ii) \\-\\ cov is the covering norm (see definition 4 in 
section 2). 

Proof. Inequalities (jHJ) follow from (JI2"]) . ( l4"3l and the upper bound 

^ p ,m s , a ^ w illf H T ifx...xlx...x5 i v IU ( 4 5) 

^x-'XlX'-xSiv'^" X 

t 

71 

constituting relation (1B18j) of lemma 5 in appendix B. ■ 

Propositions 3, 4 imply the following general statements on an Si x • ■ ■ x ^-setting LHV 
description of an iV-partite quantum state. 

Proposition 5 (a) Every N -partite quantum state p admits an 1 x • • • x 1 x S n x 1 x • ■ • x 1- 

setting LHV description; 

(b) If, for a state p on a Hilbert space "Hi (g> • • • (g> Hn, there exists a tensor positive source 

operator Tg^ x ,„ xlx ,„ xS for some n, then p admits the Si x • • • x S n x • • • x SN-setting LHV 

t 

71 
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description for an arbitrary number S n > 1 of settings at this n-th site; 

(c) If, for a state p on a Hilbert space T-L\ £g> • ■ ■ £g> Hn, there exists a tensor positive source 
operator Tg^ x .„ xS , then p admits the Si x • • • x S' n X • • • X S^-setting LHV description for 
an arbitrary number S' n of measurements at every n-th site. 

Proof. Since T[ x ,,, xl = p and ||p|| c <w = 1; from bound (jHJ) it follows that T^,,,^ x ... xl = 1 
for every iV-partite state p. By proposition 3, this proves statement (a). 

(p) 

If, for an iV-partite state p, there exists a tensor positive source operator T<f i ' x ... xlx ^ xSN > 

t 

n 

for some n, then from relation f|T3|) it follows that, for this source operator, the covering 
norm \\Tq v v1v v<? I \ cov = 1. In view of bound (HH), the latter implies ~ _ =1 

II 5ix---xlx---x6jvllcof V! vi f S!X--xS„x-xS N 

n 

for any number S n > 1 of settings at this n-th site. By proposition 3, this proves statement 

(b) . 

(p) ® 

Let an iV-partite state p have a tensor positive source operator Tg iX ,,, xS > 0. Then, for 
each n = 1, N, the operator on 'H^ 51 '' <8> • • -®1-L n <8> • • •®'H^ iV ' > reduced from T^ x ,,, xSn con- 
stitutes a tensor positive source operator ?5^ x ,,, xlx ... Xig for state p and, therefore, statement 

t 

n 

(c) follows from statement (b). ■ 

Statement (a) of proposition 5 agrees with proposition 2 of Ref.— on the LHV description 
of a general correlation scenario with setting S x 1 x • • • x 1 . 

Specified for a bipartite case (N = 2), statements (b), (c) of proposition 5 are consistent 
in view of note^ with theorems 1, 2 in Ref.—. 



VI. QUANTUM VIOLATIONS OF BELL-TYPE INEQUALITIES 

Consider a linear combination ( |26i) of averages (127|) . arising under a quantum Si x • • • x Sn- 
setting correlation scenario £ Py m sa and specified by a family ^s,a = {^( Sl ,..., Sjv )} of bounded 
measurable real- valued functions ip^ si s ) : Ai x • • • x Ajv — >■ M. 

By theorem 2, every quantum scenario £ Pi m sa admits an LqHV model and, by theorem 
1, the latter is equivalent to the existence of a bounded real-valued measure Ps pM re- 

(£ M ) 

turning all joint probability distributions P, s ^' of scenario £ Py m s A as the corresponding 
marginals. Therefore, for a quantum scenario £ p> m sa , a linear combination fl26l) of averages 
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takes the form 



E (^..^(Ai,...^)^ (46) 

Sl,...,Sjv ' 

= / E ^,... > . w) (aS' i) ,... ) A^) A* w (dA?> x ••• x dA<*> 

S1,...,SJV 

x ••• x dA$ x ••• x dA^). 



Substituting the Jordan decomposition (see remark 2) of measure u into fj4"6T) and 

£ P,M SA 

taking into the account (ITU1) . (1201) . we derive: 



inf 



£ P. M S,A 



7SUp _ yqinf "I 



^ E (V'( sl ,... lSjv )(Ai,-,A i v) 



< R sup + 
- °*s.a ^ 



P. M S,A 



P,M S A 



- 1 

(•«/• /■ jjsup einf 



(47) 



/»sup _ nini \ 
^°*S,A °*S.A^' 



where -B^J A , B^ s A are the LHV constants f )33|) and 



> 1 is the total variation 



norm of measure p e 



P, M S A 



Since inequalities (H7|) hold for each measure /i^ m ^ a returning all joint probability dis- 



tributions P f „ P ' M f\ of scenario S m,, as the corresponding marginals, we have: 



»inf _ 



< 



'S,A 

S1,...,SJV 



/ K?sup _ »inf 

(VW-,s w )(Ai>--->Ajv) ) 



P. M S,A 



7. - 1 

■sup IC P. M S,A /■ tosup _ loinf 

°*S,A °*S,A / ' 



< s sup + 



(48) 



where „ = inf u 



P> M S,A 



P. M S,A 



> 1 is the scenario parameter fT4T]) . 



Maximizing (T48]) over all possible scenarios £ Pi m sa , performed on a quantum state p and 



(M) 

Six---xSjv 

is the state parameter (T4"2"l) . for an A^-partite quantum state p and a function collection 



with outcomes in a set A, and taking into the account that sup MsA ls pMsA — Y 
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^s,a — {i J ( Sl ,...,s )}> we derive the following analogs 



B 



inf 



t(p, a ) _ 1 
x Six-xS N 1 



< (3 SU P U m 



inf \ 

S,A' 



(49) 



si,...,s n 



< R sup 4- 



T 



(M) 

5i x — xSn 



- 1 



P> M S,A 



/»sup _ »inf 
2 ^°*s,A °*S,A' 



of the LHV constraints fl32|) . Since inequalities ( I49p are non-trivial only for those ^s,a that 
correspond via (132 j) to Bell-type inequalities, we refer to (I49p as i/ie analogs of Bell-type 
inequalities for an N -partite quantum state p. 
From (SSI) it follows 



X] (^ (si ,...,, JV )(Ai,...,A iV ) 



P,M S A 



< T 



(P,A) 

Six — xSn 



By 



S.A' 



(50) 



where £>* SA is the LHV constant ( 155]) . 



Remark 3 ITie quantum constraints fffi ), ( f50j) are equivalent iff B^ 3 = —B S ^ A . For an 
arbitrary function collection ^s,A: ( fT^P =>■ ( f^Djj 6u£ no£ mce versa. In order to see a difference 
between these two types of quantum constraints for an arbitrary ^?s,a> ^ us take ^S,A f or 
which B b ~ p = 0. In this case, the left-hand and the right-hand sides of ( [T^P are equal to 



T 



(M) 



I B inf I 



T 



(M) 

St_X---xSa 



inf 



\B 



respectively, whereas the right hand side of / T50)) is given by Y^'^Lxs J 



(51) 



./Vote that, 



specifically for bipartite Bell-type inequalities with B s ~ 
two-qubit states have been analyzed numerically in Ref. 



sup 

S,A 
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0, the maximal violations by 



The following statement (proved in appendix C) shows that the quantum constraints 
( 149|) . ( 150|) are tight in the sense that the state parameter T^'^LxS represents the maximal 
violation by state p of all Bell-type inequalities (either on correlation functions or on joint 
probabilities) for a given outcome set A and settings Li x ■ ■ • x L N with Li < Si, L N < Sjy. 



Lemma 3 In / T50j) . parameter ^^ix^-xSn = su Pm sa 7£ m ^ s otherwise expressed by 



P> M S,A 



T (p,A) 

1 Six-xS N 



sup 

Ms.a,*s 



By 



S,A 



(^(s u ...,s n )(^-^n) 



P, M S A 



(52) 



V S,A' 
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where supremum is taken over all non-trivial (<B* SA ^ 0) families ^s,a — {' l l J (s 1 s )} °f 
bounded measurable real-valued functions on A = A x x • ■ ■ x K N and over all possible families 
M^a = {m1 s ™' ) } of POV measures on spaces (A n , J^An)- 



(p) 



From lemma 3 it follows that the state parameter Tg 



x-xSiv 



sup A T 



Six---xSjv> 



intro- 



duced by relation (H3"|) and discussed in proposition 3, is otherwise expressed by 



1 SiX-xSn 



sup 

A, M SiA ,* Si a, 



Sl,...,8, 



(^ 1 ,...,s JV )(Ai,-,A JV )) 



(53) 



and, therefore, represents the maximal violation by state p of all Bell-type inequalities on 
correlation functions and joint probabilities for settings up to setting S± x • • • x Sn and an 
arbitrary outcome set A. 



Proposition 6 For a state p on a Hilbert space Hi <S> • ■ ■ <8> Hn, the following statements 
are mutually equivalent: 

(a) State p admits the Si x • • ■ x Siy-setting LEV description; 

(b) Parameter T$ x „. xSjl = 1; 

(c) State p does not violate any Bell-type inequality with settings L\ < Si, ...,L N < S^ and 
outcomes in an arbitrary set A. 



Proof. Equivalence (a)-^(b) follows from proposition 3. Implication (a) =>- (c) follows from 
definition 6 of a Bell-type inequality. Let (c) hold. Then from (15B1 it follows T^ x ... xg = 1, 
so that (c) =>• (b). Thus, we have proved (a) <^ (b), (a) =>■ (c), (c) =>- (b). These implications 
prove the mutual equivalence of statements (a), (b), (c). ■ 

The following theorem introduces a general analytical upper bound on the maximal vio- 
lation T<£ x ,,, x5 by state p of all Si x • • • x S^-setting Bell-type inequalities - the maximal 
Si x • • ■ x S^-setting Bell violation for state p, for short. 



Theorem 3 For every quantum state p on a Hilbert space Hi £§> • • • <g> Hn and arbitrary 
positive integers Si, Sn > 1, the maximal Si x • ■ • x SN-setting Bell violation T<£ x ... xS > 
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1 is upper bounded by 



T (p) < inf NT (p) II (^A\ 

L Six-xS N — , . 1111 II j 5ix-x1x-x5jv W cov 



n(p) 

L Six---xlx---x 5jy 
t 



V/i 



<T inf I l'T-'(P) II 

uu l| - t 5ix---xlx---x5 JV Ml' 



T (p) Vn 
- t Si x---xlx--xS„' 



where infimum is taken over all source operators Tg^ x ,.. xlx ... xS with only one setting atn-th 

t 

n 

site and over all n = 1, ...,N and \\-\\ cov , \\-\\i mean the covering norm and the trace norm, 
respectively. 



Proof. The statement follows from relation f l53|) . proposition 4 and bound (112)) . ■ 

A. Numerical estimates 

In this section, via the analytical upper bound ( 154"]) we estimate the maximal Si x • • • x Sn- 
setting Bell violation T^ x ... xg in terms of numerical characteristics of quantum correlation 
scenarios such as a number N > 2 of sites, a number S n > 1 of measurements and the 
Hilbert space dimension c? n := dim"H n at each n-th of N sites. 

Let us first evaluate ^Si X ... x s N f° r some concrete quantum states generally used in quan- 
tum information processing. 

For the two-qubit singlet ip sing i e t = ^75 ( e i ® e 2 — ^2 ® ei), the analytical upper bound ( 154]) 
and relation (1A23)) imply 

Note that, due to Tsirelson's bound 5 and the analysis of Fine^, violation by a bipartite state 
p of an arbitrary 2 x 2-setting Bell-type inequality (either on correlation functions or on 
joint probabilities) for two settings and two outcomes per site cannot exceed \/2 - in our 
notation T^^ 1 '^ ^ < \/2. The maximal violation by the two-qubit singlet of all correlation 
Bell- type inequalities is given^ by the Grothendieck's constant Kq(3) of order 3 and it is 
known!^ that V2 < K G (3) < 1.5164. 

Consider also the maximal Bell violations for the N-qudit Greenberger- Home - Zeilinger 
(GHZ) state 

Vd j=1 
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and the generalized iV-qubit GHZ state 



(gen) 



simp \1)® N + cos <p |2) 



2 ) 



(57) 



where j = 1, d, are mutually orthogonal unit vectors in C d . For each of these states, 
the trace norm of the source operator (1A20I) is upper bounded by 

-{p d ) 



IX S2 X ••• XSjy 

r(p^ 9en) ) 

IX S2 X ••• X5jy 



< l + 2 7V - 1 (d-l 



< 1 + 2 W 1 I sin 9? cos 9? I 



and, in view of relations (1351) . (IA90 . the analytical upper bound (154"I) implies 

< min{(25 - l)^" 1 , 1 + 2 N ~ 1 (d - 1)} 



L Sx---xS 



and 



T 



(4 9en) ) 
5x-x 



< 1 + [min{S N -\d}-l\ 



s < min {(2S - l) N -\ 1 + 2 N ~ X |sin2^|} 



(58) 



(59) 



(60) 



< 1 + 2 N ~ l \sm2p\ • 

The second lines in ([59]), flSOD are due to relation (25 - l) 7 ^ 1 < 2 N ~ 1 (S N - 1 - 1) + 1 that 
can be easily proved by induction. 

For an arbitrary Af-partite state p, the general analytical upper bound (1541) implies the 
following new numerical upper estimate. 

Theorem 4 For every state p on a Hilbert space Hi <g) • • • <g> and arbitrary positive 
integers S%, Sn > 1, the maximal Si x • • - x SN-setting Bell violation satisfies relation 



T^ x ... x5jv < min{^, 9 N , } 



(61) 



< 1 + 2 Af " 1 
where d n = dim"H n , n = 1, N, and 

- 1 



mm 



di ■ . . . ■ dj\[ Si ■ . . . ■ S 
maXy, dn max^ S r 



N 



^ 1 + 2 *-i (di-...-d N 







N 



-1 



\JV-1 



N-2 

min V(-l) fc 2 7V " 1 - 

{ni....,njv-i} ' 
C{1,...,JV} fc=0 ^i^-^jv-l-fc' 

UjG{n 1 ,...,n N _ 1 } 



(62) 



E 



5, 



'.'1 
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If, in particular, d n = d, S n = S, Vn, then 



T 



(p) 

S x ... x 



s < min {(25 - I)*" 1 , 2 N -\d N - 1 - 1) + 1} 



(63) 



< 1 + 2 



N-l 



(min{S, d})" -1 - 1 



Proof. From bound (|54p it follows 



y(p) <-- 
1 S 1 x-xS N — 



inf 



J S 1 x-xlx--xS JV > vrt 

t 



< inf 

1 S 1 x---xlx---xS N ^ n 



<min|||r^ . „ |U„ \\r ip) 



where 

5] X-Xl X---X 5jy 



and 

.S-j X---X 1 X---X Sjy 

t 



t 

n 

\\T (p) Hi 

11 S 1 x---xlx---xS N M± 

n = l,...,iV}, 
are the specific source operators constructed for an 



arbitrary iV-partite state p in appendix A. Taking into the account the upper bounds (IA9I) . 
( IA2ip for these source operators and also relation (2S — l)^ -1 < 2 N ~ 1 S N ~ 1 — 2 N ~ 1 + 1, we 
come to relation f lBTj) . implying, in turn, f l6"3"j) . ■ 

Note that estimate flBTj) implies that T^|,„ xlxS xlx ... xl = 1 for every Af-partite state p. In 
view of statement (b) of proposition 3, this result agrees with statement (a) of proposition 
5. 

Theorem 4 implies. 



Corollary 2 (a) For an arbitrary N -partite quantum state, violation of a Bell-type inequal- 
ity (either on correlation functions or on joint probabilities) for S settings per site cannot 
exceed (2S — l) 7 ^ 1 even in case of an infinite dimensional state and infinitely many out- 
comes. 

(b) For an arbitrary state p on (C. d )® N , violation of a Bell-type inequality (either on correla- 
tion functions or on joint probabilities) is upper bounded by 2 N ~ 1 (d N ~ 1 — 1) + 1 independently 
on a number of settings and a number of outcomes at each site. 



Let us now specify the general A^-partite upper estimate fl6"2"|) for N = 2,3. 
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Corollary 3 For every bipartite state p and arbitrary positive integers Si, S% > 1, 

T^ x52 < 2wm{Si,S 2 ,d 1 ,d 2 } - 1. (65) 
If, in particular, d\ — d 2 = d, Si = S 2 = S, then 

T i s p) xS <2mm{S,d}-l. (66) 
For every tripartite state p and arbitrary positive integers Si, S 2 , S3 > 1, 

^Ux5 3 < min ( min (4S ni S n2 -2(S ni +S n2 ) + l), 4-^^-3) (67) 
U»i,n 2 }c{i,2,3} max n cs n j 

f S1S2S3 did 2 d 3 \ 

< 4 mm < — , — V - 3. 

[ max n S n max n rf n J 

//, in particular, d n = d, S n = S, Vn, i/ien 

rg 5x5 <min{(2S-l) 2 , 4rf 2 - 3} (68) 
< 4(min{S, d}) 2 -3. 

We stress that, in contrast to the bipartite and tripartite numerical estimates found in 
Refs.— >22'~— up to unknown universal constants, our bipartite and tripartite numerical upper 
estimates fl65|) - (16"8"]) are exact. 

B. Discussion 

For bipartite and tripartite correlation scenarios with a finite number of outcomes at each 
site, the numerical estimates on the maximal Bell violations have been recently presented 
in Refs.— ~— . The results of corollary 3 indicate. 

• Our exact bipartite upper estimate ( l6"5"j) improves the approximate bipartite estimate 
^ min{<i, S} found in Ref.— (theorem 6.8) up to an unknown universal constant. 

• The bipartite upper estimates (in our notation) 

T { S P ; X %<2K G + 1, if 14 = 181 = 2, (69) 

T^ A j 2 < 2 |^|i3|(^ G + !)-!, VL4|,|B|, 
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derived in theorem 22 of Ref.— for arbitrary S\,S2 > 1 and numbers |*4.] , \B\ of 
outcomes at Alice's and Bob's sites, are improved by our bipartite upper estimate 
(J65D if 

min{di, d 2 } < K G + 1, (70) 
min{di,d 2 }< \A\ \B\ (K G + 1), 

respectively. Here, Kq is the Grothendieck constant, and it is known^^ - that Kq G 
[1.676,1.783]. 

• From our exact tripartite upper estimates (1671) . (168]) it follows that violation by a 
tripartite quantum state of a Bell-type inequality for S settings per site cannot exceed 
(2S — l) 2 . Therefore, the tripartite lower estimate >z yfd, found in theorem 1 of Ref.— 
for violation of some correlation Bell-type inequality by some tripartite state on C d £g> 
C D (g> C D , is meaningful if only a number S of settings per site needed for such a 
violation in the corresponding Bell-type inequality obeys relation 

(2S - l) 2 h Vd. (71) 

Thus, for an arbitrarily large tripartite violation argued in Ref.— to be reached, not 
only a Hilbert space dimension d but also a number S of settings per site in the 
corresponding tripartite Bell-type inequality must be large and the required growth of 
S with respect to d is given by (ITTT) . 

VII. CONCLUSIONS 

In the present paper, for the probabilistic description of a general correlation scenario, we 
have introduced (definition 5) a new simulation model, a local quasi hidden variable (LqHV) 
model, where locality and the measure-theoretic structure inherent to an LHV model are 
preserved but positivity of a simulation measure is dropped. 

We have specified (theorem 1) a necessary and sufficient condition for LqHV modelling 
and, based on this, proved (theorem 2) that every quantum correlation scenario admits an 
LqHV simulation. 

Via the LqHV approach, we have constructed analogs (Eq. (1191) ) of Bell-type inequalities 
for an iV-partite quantum state and found (theorem 3) a new analytical upper bound on the 
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maximal violation by an iV-partite state of all Si X • • • x S^r-setting Bell-type inequalities 
- either on correlation functions or on joint probabilities and for outcomes of an arbitrary 
spectral type, discrete or continuous. 

This analytical upper bound is based on the new state dilation characteristics (definitions 
1-4 ) introduced in the present paper and this allows us: 

• to trace (propositions 5, 6 ) iV-partite quantum states admitting an S\ x • • • x Sn- 
setting LHV description; 

• to find the exact numerical upper estimates (Eqs. fl55|) . fl59l) . fl60|) ) on the maximal 
Si X • • • x S^r-setting Bell violations for some concrete iV-partite quantum states used 
in quantum information processing; 

• to find (theorem 4 ) the exact numerical upper estimate on the maximal Si x ■ ■ • x Sn- 
setting Bell violation for an arbitrary iV-partite quantum state, in particular, to show 
(corollary 2) that violation by an iV-partite quantum state of a Bell-type inequality 
(either on correlation functions or on joint probabilities) for S settings per site is upper 
bounded by (2S — l)^ -1 even in case of an infinite dimensional quantum state and 
infinitely many outcomes. 

Specified (corollary 3) for N = 2,3, our exact iV-partite numerical upper estimate (Eq. 
fl6T|) ) improves the bipartite numerical upper estimates in Refs.— ^ and clarifies the range 
of applicability of the approximate tripartite numerical lower estimate in Ref.— . 

Acknowledgement 1 / am grateful to the organizers of the Operator Structures in Quan- 
tum Information Workshop, held at the Fields Institute in Toronto on July 6 -10, 2009, 
where part of this paper covering a bipartite case was presented. I am also thankful to Pro- 
fessor A. M. Chebotarev for the useful discussion. 

Appendix A: proofs for section 2 

In this appendix, we prove proposition 1 on the existence of source operators for an N- 
partite quantum state and introduce some source operators which are different from those 
constructed in the proof of proposition 1 and are needed for our consideration in section 2 
of appendix B. We also prove lemmal on the properties of the covering norm. 
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1. Proof of proposition 1 



For a state p on a Hilbert space H 



1-Lni consider its decomposition 



Iromi fcfci = ^ a «^m-^mi,..fa' «j > 0, ^^a» = l, 



(Al) 



(0 

p>m...fc 



m,...,k 



in orthonormal bases {e™ G n — 1, ...,iV. 

Let X = 2. For a bipartite state p, denote by p n the reduced state on "H n , n = 1, 2, and 
introduce on a Hilbert space ftf 5l <g> ftf 52 the self-adjoint operator 



T 



(p) 

SixS 2 



Vmmi,kki 



sym 



(A2) 



sym 



(S 2 - l)pf 51 ® af 52 - (Si - l)of 51 ^ - 052 



P2 



-(S 1 -l)(S 2 -i)<Tf 3l ®af 8 ', 



where a n is a state on H n and notation [-] sym means symmetrization on l-i® Sn . For example, 
[Xi ® X 2 ]^ m := Xi <g> X 2 + X 2 (g> X±. It is easy to verify that (1A2j) represents an Si x SV 
setting source operator for state (lAlj) specified with X = 2 and the trace norm of this source 
operator satisfies relation 



1 < 



< 2S 1 S 2 - 1. 



(A3) 



Let X = 3. For a tripartite state p, in addition to the above notation p n for the reduced 
state on H n , denote by p ln the reduced state on %i ®% n . For short of notation, we further 
take one of settings to be equal only to 1, say S\ — 1. Introduce on a Hilbert space Hi eg) 



Kf 52 ® "Hf 53 the self-adjoint operator 



T 



(p) 

lxS 2 xS 3 



Immi ,.i.hMi I e m^ ) ( e m! I ® l e j ^)( e ji^l ® a 



5(Sa-l) 



sym 



(A4) 



4 3, )<4?l 



a. 



5(53-1) 



sym 



®S3 

3 



(Pl,3) 
r lxS 3 



- (S 2 - IV ® erf 2 ® / - (S 2 - 1)(S 3 - l)pf Sl ® 



<7 



®5 3 
3 ' 



where: (i) r[ P ^ is the 1 x S n -source operator (IA2j) specified for the reduced state p 1 n on 
"Hi ®T-i n ] (ii) notation in the third line means operator derived by insertion of term ®0"f 52 ® 



into each term of the tensor product decomposition of the source operator r^sf for state 
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p x 3 on Hi <g> "H 3 . It is easy to verify that (1A4j) represents an 1 x S 2 x ^-setting source 
operator for state flAip specified with iV = 3. Substituting ( 1A2j) into (1A4[) . taking into the 
account (jSJ) and evaluating the negative part of the self-adjoint operator (IA4I) . we derive 



T 



IXS2XS3 



< 1 + 2(S 2 - l)S 3 + 2(S 3 - 1)S 2 
= AS 2 S 3 - 2{S 2 + S 3 ) + 1 

< 4S 2 S 3 - 3. 



(A5) 



Let N = 4. For a quadripartite state p, denote by p n , p nin2 , P nin2n3 ^ ne reduced states on 
Hn, U ni ®Un 2 and U ni ®H n2 ®U nA , respectively and consider on Ui ®U 2 S2 ®Uf S3 ®nf Si 
the self-adjoint operator 



T 



(P) 

1 x $2 x 53 x 54 



,S y^ J T lm^ jji&hkki l e m^)( e mi I ® l e j ^)( e ji^l ® a 

[h (3, >(i 3, i « 



sym 



5(^3-1)' 



J sym 



4 4) ><4>*i 



5(^4-1)' 



J sym 



(Pl,3,4) 
r l x S 3 x 5 4 



(S 4 - l)rj 



(Pl,2,3) 
xS 2 xS 3 



a 



®S4 



- (S 2 - 1)/® <7 



®S 2 
2 



/ ® / 



(Pl,2,4) 
r lxS 2 xS 4 



- (5a- 1)/® /®af 53 ® / 
-(5 3 -l)(5 4 -l)rM® 



<2>S 3 ^.®5 4 



£7 o 

(Pl,4> 
T lxS 4 



- (S 2 - 1)(S 3 - iy ® a® 52 ® af 53 (8) / 



Ol,3> 
r lxS 3 



(5 2 - 1)(5 4 - 1)/ <8> af S2 g> / ®a 



- (5 2 - 1)(S 3 - 1)(S 4 - l)pf Sl ® af 52 ® a 



4 ? 



(A6) 



where r^sf is the 1 x ^-setting source operator ( 1A2I) specified for the reduced state p l n 
and rlxs^'xln is the source operator (IA5|) for the reduced state Pi nim - It is easy to verify 
that 0A6p represents an 1 x S 2 x S 3 x ^-setting source operator for state ( \Al\i specified 
with N = 4. Substituting flA2j) . ()A4j) into ( 1A6j) and evaluating the negative part of the 
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self-adjoint operator (1A6[) . we derive 



-to 

1 x S 2 x S 3 x S 4 



< 1 + 2(S 2 - 1)5 3 S 4 + 2(S 3 - 1)S 2 S 4 

+ 2(S 4 - 1)S 2 S 3 + 2(S 2 - 1){S 3 - 1)(S, - 1) 
= 8S2S3S4 — 4(^5*3 + S 2 S4 + 
+ 2{S 2 + S 3 + S 4 )-1. 



(A7) 



Thus, in view of (JMj), (TA"S]h (jATj) . 



-to 

lxS 2 



-to 

lxS 2 xS 4 



-to 

1 X S2 X S3 X S4 



< 25 2 - 1, 

< 4S 2 S 3 - 2(5 2 + S 3 ) + 1 < 4S 2 S 3 - 3, 

< 8S2S3S4 — 4(S 2 S 3 + ^5*4 + S3S4) 
+ 2(S 2 + S 3 + S 4 ) - 1 < 8S 2 S 3 S 4 - 7. 



(A8) 



All these bounds are tight in the sense they imply \\t[ p J xi || = 1. The generalization to an 
A-partite case with setting S% = 1 is straightforward and gives 

JV-2 

^^(-1)^2^ S^-.-.-Sn^ (A9) 



-to 

lxS 2 x ■■■ x5jy 



fc=0 



ni^...^n ]v _i_fc, 
rij=2,...,N 



+ (-1 



,N-1 



-to 



lxSx-xS 



< (25-1 



,iV-l 



The constructed source operators (1A2I) . (1A4I) . (1A6I) . (1A9I) prove the statement of propo- 
sition 1. 



a. Other examples of source operators 

For our consideration in section 2 of appendix B, let us also construct source operators 
of a type r different from type r in Eqs. fTA~2"]) . fTA"4~|) . fTA~6]) . (jMJ)- Denote d„ := dimT^, 
n = 1 , . . . , A, and assume that max n <i n = d\. 

Let N = 2 and i[) & Hi (g> "H 2 . For a pure bipartite state |^)(^|, consider its Schmidt 
decomposition 

mm = ® isf x^fi, ^ > o, ^ 2 = 1, (aio) 
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where the sum is taken over = 1, d 2 and {g, } is an orthonormal base in T-L nj n — 1, 2 



Introduce the self-adjoint operator 



lxS 2 



(2,S 2 ) 



(All) 



where 



2W. 



(2,5a) 



/, (2) . (2) w (2) . (2)|\® 52 / 1 (2) (2) w (2) (2) A 



®5 2 



jji y 



2 5 2 

(|jf + + is. 



(2) 



2 5 2 

'\ (2) • (2)w (2) • (2) I x 052 



(A12) 



2 s 2 



2 s 2 



(2,S 2 ) 



It is easy to verify that flAllj) represents an 1 x ^-setting source operator for state (lAlOj) 
and 

2 



lxS 2 



1 < 1 + 2 E^ 1 = 2 (S^) - 1 



(A13) 



<2do-l 



for any S2 > 1. By convexity, for an arbitrary bipartite state p = J2i Q; j|V'i)(V'J; operator 



T 



(p) 
lxS 2 



(I^X^I) 

lxS 2 



represents an 1 x ^-setting source operator and 



dp) 

lxS 2 



< 2d 9 - 1, V& > 1. 



Let N = 3. For state (IA1I) with N = 3 introduce the self-adjoint operator 



r(p) 
1 x S 2 x S 3 



E 



Vmmi ,jjx...kki I 



(A14) 



(A15) 



where W, 



» 



<8>>S„ 



r(n,S„) 



and operator is defined by (IA12j) via replacements 

2 — > n, 5*2 — > S n , — > . It is easy to verify that (|A15[) is an 1 x S 2 x ^-setting source 
operator for state flAlj) with AT = 3. Splitting (1A15|) into four sums 



-(p) 



3 (2)\/X2)|N(g)5 2 



3 ( 3 )\/J 3 )h«>53 



(A16) 



,(2)\/ (2)h®Si 



W, 



(3,S 3 ) 



hj^juk^kx 



r(2,5 2 ) 



r(3,S 3 ) 
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where, for any index i, 



— Vc (,) e {1) \\6 (t) 



Pi' 



E 



1/2 



(**) 1 2 



mj'fc I 



W\2 
jfcJ 



and taking into the account that 



(n,S„) 



1 W (n,5 ' l) 



<2, 71 = 2,3, 



we derive 



lxS 2 xS 3 



< 4cMs - 3 



(A17) 



(A18) 



(A19) 



for any S 2 , £3 > 1. 

The generalization of flA15j) . flA19[) for A > 4 is straightforward and gives the source 
operator 



IXS2 X ■••XSjy 

with the trace norm 



/ ^ /mmijj'i roil tn /\ mi I ^ jji 



w, 



(N,S N ) 



kki ' 



(A20) 



r 



(/>) 

1 X S2 x "~ X $N 



<2 N -\d 2 - ...-d N -l) + l. 



(A21) 



6. Source operator for the singlet 



For the two-qubit singlet ip s i ng i et = ^75 ( e i ® e2 — e 2 <S> ei), the self-adjoint operator 



T'iPsinqlet 
1 



ix2 } = — |ei> <ei| ® |e 2 )(e 2 | ® |e 2 )(e 2 | + -|e 2 )(e 2 | ® | ei )<ei| <g> |ei><ei| (A22) 

- ||ei><e 2 | ® |e 2 )( ei | ® - \\e 2 ){e l \ ® \e 1 ){e 2 \ ® 

- ||ei><e 2 | <8) -y <g> |e 2 )( ei | - ^|e 2 )( ei | ® ® | ei )(e 2 | 

on (C 2 ) 03 represents an 1 x 2-setting source operator. Calculating the eigenvalues of this 
source operator, we derive Ai )2i3i4 = 0, A 5j6 = and A 7i8 = Hence, 



{.P singlet ) 



1x2 



Vs. 



(A23) 
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2. Proof of lemma 1 



Property 1. The first left-hand side inequality of property (1) is trivial. The last right- 
hand side inequality is due to the fact that, for each W G 7gi®...®g , operator \W\ is one of 
its coverings, so that, by ffTTj) and definition of the trace norm, we have 



\W\\ cov <tr[\W\) := \\W\l, 



In order to prove 



sup 

X S =X}, ||X,||=l,Vj 



tr[W{Xi®"-®*m}]|< \\W\\ cov , 



(A24) 



(A25) 



where supremum is taken over all self-adjoint bounded linear operators Xj, j = 1, N, with 
the operator norm \\Xj\\ = 1, let us represent operator W G Tg \®...®g y i a decomposition 
([7]) with an arbitrary a trace class covering W cov . We have 



sup | tr[W{Xi 



® X m }} | < - sup | tv[(W C0V + W){X 1 <g> • • - ® X m }] | 



(A26) 



+ -sup| tr[(^-^){Xx 



X m }] | . 



Applying to each term in f lA26j) the spectral theorem 



Xi. 



AE Xfc (dA), 



(A27) 



where K Xk is the spectral (projection- valued) measure on (R, £>k) for each self-adjoint 
bounded linear operator X^, k = 1, ...,m, and taking into the account that W cov ±W > 
and spectrum a(Xk) C [—1,1] for each fc, we derive 



sup | tr[(W cov ±W){Xi®---® A m }] | 



(A28) 



< 



sup 



y Ai ..... A m tr [(W„» ± WO {E Xl (dA x ) <8> • • • ® E Xm (dA m )}] 



< tr [Fy cot , ± W] 



Substituting (IA28[) into (IA26I) . we have 



sup |tr[W{Xi 

x s =x;, \\\x 3 \\=i,vj 



X m }] | < tr [W a 



(A29) 



for each trace class covering W cov of an operator W G 7g^„,^g ■ This relation and definition 
ffTTj) of the covering norm imply inequality (IA25|) . 
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Property 2. If a self-adjoint trace class operator W is tensor positive, then tr[W] > 
and W is itself one of its possible coverings. Therefore, by ( ITT]) . ||W^|| col , < tr[W]. This and 
property (1) of lemma 1 prove property (2). 

Property 3. Let W re d G "Tg S k %>...®g k .i 1 < ki < ■ ■ ■ < kj < m, j < m, be the self-adjoint 
trace class operator reduced from a self-adjoint trace class operator W on a Hilbert space 
Gi <8> • • • ® £? m - The left-hand side inequality in property (3) follows from property (1) and 
relation tr[W re d] = tr[W]. 

Further, if (Wcot,)^ G Tg kl ®-®Qk- ls the operator reduced from a trace class covering 
PVcot, of operator G 7gi®-®g > then (W C(W ) r e<i is one of trace class coverings of operator 
Wred £ 'l^h%---®G k .- Therefore, for each W G Tg*®...®g m , we have the following inclusion 

{W cov ) red G 7^ 



red) cov ^ Tg kl ®...®g k . 

In view of definition ( ITT]) of the covering norm and relation tr [(W / CO t))red] = tr[W ccw ], inclusion 
dA30]) implies 

||W rerf || C(W = inf ti[(W red ) cov ] (A31) 

(M / red) cot ,S7g J . i ®...®g fc 

< inf tr [(W^,) 

WcotrG7ff 1 ®...8ff m 

= inf tr[^] = ||^||_ 
for each W ret j reduced from an operator G 7g^...^g m - This proves property (3). 



Appendix B: proofs for section 5 

In this appendix we prove lemma 2 and find also some upper bounds needed for our 
presentation in section 5. 



1. Proof of lemma 2 



For short, we further omit the below indices p, M S) \ at notation £ Pi m sa an d denote by 
M.£ := {p} the set of all normalized bounded real- valued measures p, each returning all 
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distributions P, si . of scenario S as the corresponding marginals and denned on the direct 
product space (Af 1 X • • • X A S N N , Tff 1 <g> • • • ® J 7 ^) := (A', J 7 ')- 

If scenario £ admits an LHV model, then, by corollary 1, in set Ais, there is a probability 
measure. Since the total variation norm of any probability measure is equal to 1, parameter 

IB = hlf^eA/jg ||/i||TOir = 1- 

Conversely, let j s = 1. As it is shown in section 2 of this appendix, for a quantum 
scenario, set contains more than one element, so that, in view of convexity of Aig, this 
set is infinite. From relation 

7g = inf{|H„ ar \fteM s } = l (Bl) 

it follows that, in set Aie, there is a sequence {/x m }, for which ||/vll«ir 1 as m -)• oo and 
which is bounded in norm || • \ \ var . Note that, equipped with the total variation norm ||-|| mr , 
the linear space 5(A',.F') of all bounded real-valued measures on the measurable space (A', J 7 ') 
is Banach. Therefore, there exists a subsequence {/ifc m } Q {^ m } an d a measure Ji G 3ya',.f') 
such that 

lim / /(A')^(dA') = / fWZ(d\% \\Jl\lar < 1, (B2) 
m->oo y y 

for all Borel measurable bounded real-valued functions / on (A', J 7 '). Denote by Pr slj ___ jSN ^(dXi x 

• • • x dAjiy ) the corresponding marginal of measure ju(dA^ x • • • x dAp 1 x • • • x dA^ x • • • x 

dA^). Specifying flB2) with the indicator function x F ((Af l} , A& w) )), F G J 7 ! ® ■ ■ • <g> T N , 

and taking into the account that fi km G A4g, Vfc m , we have 

P(Z..,sJ F ) = P( Sl ,...,s N) (F) (B3) 

for all sets F 6 Ji ® • ■ • ® Fn and all tuples (si, sjy). Thus, P( S1) ..., SJV ) = PA s % for all 
joint measurements (si, sjv) of scenario £. The latter means that the bounded real- valued 
measure /x is normalized and returns all probability distributions P) s of scenario £ as 
the corresponding marginals. Hence, // belongs to set A^£-. In view of relation (|2"U|) and 
the second relation in (IB2I) . the normalized bounded real- valued measure Jl G .Mg is a 
probability measure. By corollary 1, this proves lemma 2. 

2. Some upper bounds 

For the evaluation in proposition 4 of the state parameters (1421 . (I4"3"j) . in addition to 
measure (|40p . let us also consider another example of possible measures /% pM in (|4ip . 
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Let ^ix5 2 x---x5jv be an 1 x 5*2 x •■■ x S^-setting source operator for state ponHi® - • •®%at 
(see definition 1 in section 2). introduce the following collection of normalized bounded real- 
valued measures 



K[T^ S2X ... xSn {M^\d\[^) ® M«(dA«) ® • • ■ ® Mf)(dAf 2) ^ 



(B4) 



® - • • ® Mff(dA^) ® ■ • • ® M^CdA^)}], 



si 



1, ^i, 



where each si-th measure returns as the corresponding marginals all joint distributions 
Pr 8 * MS grf) °f scenario ^ p ,m sa with measurement si at site n = 1. 



For an arbitrary trace class covering (T{£ s , xS ) cov of a source operator 7 T ^ 5 S v 
(see definition 3 in section 2), decomposition (|7j) implies the following representation 



1 

2 tr 



r+{ 



(«i), 



M«(.) 



} 



1 

2 tr 



r - { Mi' 



(«)/ 



M«(-) 



} 



of each si-th measure (lB4p via two positive real-valued measures, where 



T 



rp{P) 

1xS2X---X-Sn 



(p) 

xS 2 X---xSjv 



> 



(B5) 



(B6) 



are tensor positive trace class operators. 

As it is discussed in the proof of theorem 2 in Ref.— , for a positive measure v on a direct 
product space (Ai, J~Xi) x (A*2,J~a 2 ), eac h positive measure v{B\ x •), Bi E J-\ 17 on (A 2 , F\ 2 ) 
is absolutely continuous^ with respect to the marginal measure v(Ai x •) on (A2,J 7 a 2 )- 
Hence, for each of positive measures in decomposition (IB5|) . the Radon- Nykodim theorem^ 
implies representation 

r^M^dA? 15 ) ® M? } (dA^) ® • • • } 

a^(dAS Sl) I A?*, ...) tr [t±{% ® M«(dA«) g 



tr 



(B7) 



} 



via conditional probability measure a«^(- | A 2 , •••) and marginal 

tr [^{I^gM^dA?)® •••} 

1B7D it follows that the normalized bounded real-valued measure 

ASi) 

V . . . \s 



From (1B4 



(p,M SjA ) 

1 X S 2 X ■ ■ ■ X Sjy ' ^ ' 1 X S 2 X • • • X Sjy 



x • ■ ■ x dA] v ; x ■ ■ • x dAjy 



1 

2 



(dAS 1} x ■ ■ ■ x dAj' 
J] «W(dAS Sl) |A^,-)) tr[r + {l Wl ®M«(dA«) 



(B8) 



(B9) 



,si=l,..,S'i 



1 

2 



J] a( 7 )(dAS Sl) |A«, ...) tr [r" { l Hl ® M«(dA«) 



.si=l,..,Si 
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returns all distributions (13"9"|) of scenario £ p ,m sa as the corresponding marginals. 



Lemma 4 For arbitrary source operators T^ x ,^ xSn , T^ s x xS for state p, the total vari- 
ation norms of measures ( f^OP , liB9\) satisfy relations 

i< ^ SA) < rw.„ x5w (bio) 

1 S 1 x---xS N cou 

war 

and 



1 < 



inf 

r (p) 

1 x S2 x ■ ■ ■ x Sjy / 



. (P,M S , A ) 

A X t (p) / t (p) X 

IX S2 x ■■■ xS N > V J lx S2 X ■■■ X Sjy >cov 



< 



yO) 

1 X S2 X ■ ■ ■ X SjY 



(Bll) 



/or ei>en/ collection M^a of POV measures on A and an arbitrary outcome set A. Here, infi- 
mum is taken over all trace class coverings (T^ ) of a source operator 

J v 1 X S 2 X • ■ • X S N ' cov J r IX S' 

and \\-\\ cov is the covering norm (see definition 4 in section 2). 



<S? x ■■■ x5j\ 



Proof. In view of ((7|), consider the decomposition of a source operator 



T, 



(p) 



Six...xSn 



5fa 



(P) \ j_rp(p) 

S-lx-xSnJcov -r J 5 lX -x5iv 



5 fa 



(p) \ 
Six-xS N Jcov 



T, 



(p) 



Six — xSn 



(B12) 



00 



via tensor positive operators (Tcf 



(p) 



where (T, 



(p) 

Six---xSn J c 



is an arbitrary 



trace class covering of T<f^ x _ xS , see definition 3 of section 2. 
Substituting (1B12[) into (jlDjl . we represent measure p 



(p,m s ,a) 
r (p) 

[ Si x-xS. 



as the difference of two 



positive measures and, for each of these measures, the total variation (fl8l) is upper bounded 
by 



1 



-tr 



T, 



(p) 



Six---x5j\ 



(p) 



Si x-xSa 



(B13) 



Therefore, for the total variation ffTgj) of the normalized measure p 



(p,M s , A ) 

'Six— XSjv 



we have 



1 < 



,,(P,M S , A ) 

- I S 1 x-xS JV 

^(P) \ 



< tr 



<p) 



5iX---xSjv 



(B14) 



for every trace class covering (T^ > x ... xS ) cov of a source operator T^ x _ xS . Relations (jEEJ), 

(HID imply bound f lBTOj) . 

Quite similarly, for measure (IB9p . 



1 < 



(( (P,M S , A ) 



-.(p) 



1 x S2 x ■ ■ ■ x Sjy ' y 1 x 5 2 x ■ ■ ■ x S N j 



< tr 



1 x S2 x ■ ■ ■ x 5jy 



(B15) 
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for each trace class covering (T^ 

° v 1 x £ 

ffTTT) we derive 



S 2 x---xSjv 



) cov of a source operator Tj^ s ^ xSn - From (1B15I) . 



1 < 



inf 

r (p) 

1 x S2 x ■ ■ ■ x 5jy y 



.ApMsa) 



1 X S2 x '" x Sjv ' V 1 x ^2 x '"' X $N 



< 



J'(p) 

1 X 52 X ■ ■ ■ X Sjy 



(B16) 



representing relation (1B11[) . 



Generalizing measure f lB9j) to the case of a source operator T^ x _ xlx ,.. xS with 5 n = 1 



at an arbitrary n-th site, similarly to bound (IB lip , we have 



1 < 



inf 



frp(p) 

l J SiX"-xlx-xS 
t 



II (pMs,a) 11 

Wr^rp(p) ( rp(p) n Wvar 

jJcov t T 



(B17) 



< I IT, 



(p) 



Six—xlx—xSx Wcov- 
t 



Bounds (IB13I) . flBl Tj) allow us to evaluate the scenario parameter 1e pMsA defined by 
relation fl4"TT). 



Lemma 5 For each Si x • • • x SN-setting correlation scenario £ p ,m S a on a p on a 
Hilbert space ~K\ <S> ■ ■ ■ <8> Hat, 



'v^ < inf IIT ( ' P ' 1 II 

l£„M c . — / . M- L 5lX---XlX---xSjv Ileal! - 

i S 1 X-XlX-XSjyl V " A 

t 



(B18) 



Proof. From ED, flBTUj) . (IBTTjl . it follows 



7c < inf <i I IT. 

'tp,M S ,A — 



(p) 



1 1 cov ' 






± SiX-xS N 



, n= 1,...,N , . 



(p) 



where infimum is taken over all source operators Tg 



n = 1, ...,iV. For each n, the set {T<J^ x ... xlx ... xS } of all source operators Tg P i ' x ,,, xlx ... xSN for 



x---x5jv' 5iX---x1x---xSjv 
t 

n 



(B19) 
and over all 



(®Si) 



®H { £ N \ each 



state p includes the set of all source operators on H\ 
reduced from some Tg^ x ,„ xS , as a particular subset. This inclusion and relation (fill) imply 



inf 



r (p) 

t 



|7^(p) || < 

\ ± Six — xlx — xS N \\c°v — 

t 



inf 

r (p) 

l Si x-xS» 



(p) 



Six — xSn 



(B20) 



for each n. Taking this into the account in (IB 19j) . we prove (IB18I) . 
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Appendix C: proof of lemma 3 in section 6 



Constraints 



imply 



sup 

*S,A.S* S A #0 



E (Vw., v )(Ai,-,A 

U *S,A 



< 



7£. 



P,M SjA ' 



(CI) 



where supremum is taken over all non-trivial (i. e. B^ s A 7^ 0) function collections \I/,s,a- For 
short of notation, we further replace £ Pi m s ~* £s,A- in order to prove 



sup 



By 



— E (V'(. 1 ,... I . w )(Ai,-,Aj V ; 



S,A 



8l,...,S, 



£s,A 



7£ S ,A ' 



(C2) 



we note that, by introducing variables £ = /x f (f2) > 0, the parameter 7g gA , defined by 
OH]), can be otherwise expressed as 



P, 



(£s,a) 



S,A 







C+ TD (t S,A > 



(03) 







Vs n , Vn}. 



As it is specified in section 3, we consider only standard measurable spaces. In this case, 
(A n , J"a„) is Borel isomorphic to some measurable space (X n , Bx n ), where X n E £>r is a Borel 
subset of K and Bx n '■= B^ D X n is the trace on X n of the Borel a- algebra £>s on R. We have 
two major cases. 

(a) Discrete case. Let, for a correlation scenario £s,a, each outcome set be finite: A n = 
{A^ fe ™\ fc n = 1, K n < 00}. Then 



Sl, 



(C4) 



k\,...,k N 



where (3 



(fci,...,fcjv) 



are real numbers. Hence, in ( 101)) . supremum over ^/s ; a reduces to supremum 



over families )} of real numbers and equality (102)) follows from ( 101)) and (I03P by 

the linear programming (LP) duality. This proof is similar to the proof of theorem 17 in 
Ref.— for a bipartite case with a finite number of outcomes at each site. 

Let now, for a correlation scenario £s,a, every outcome set A n be inifinite but countable: 
A n = {X^ n \k n = 1,..., K n , ...}. Consider collections ^gf K \ of bounded measurable real- 
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valued functions 



<*2„)( A i"- A *) : = E < 1 ;:3^) } ( A i)-----x { ^) } (Aiv), (cs) 



fci,...,fcjv 



specified via tuples K := (ifi, iCjy) of positive integers Ki,...,Kn < oo and families 
B is := {B(° u '"' 1 ° N } G M, s„ = l,...,S n , k n = l,...,K n \ of real numbers. For each of these 
function collections, the expression 



E «'*,(A..-.A W ) 



S1,...,S. 



£s,A 



(C6) 



? (fcl,...,fcjy) o(^S,a) 



{A? l} } x ... x {A^ ; } 



sij.-.jS^y, ki,...,k N 

is similar by its form to (1C4j) . so that relation (IClj) and the proof in the above case imply 



7£ 



, A > sup 



> sup 

0K, K 



— E (^t!l,(^.-.v 



B,, ; . Is , 



S, A *1> — l a jv 



N, 



£s,A 



(C7) 



(K) 

sup7 fs , A , 

A' 



where, for each K = (Ki, Kn), 



ifsl ■= inf { e + + T I ^ > 0, e + - T = 1, 3£JS\ ^'a : (C8) 
P ^,s n) (RJ x ■ ■ • x KJ) = ({A<*>} x ... x {A<**>}) 

x • • • x {A^}) , fc n = 1, K n , s n = 1, S n }. 



(c lhv\ 
t~ p { S,A ) 

S r ( Sl ,...,s n ) 



From (IU31), (jCS|> it follows that 7 gp < 7^ < 7e if < L X ,...,K N < L N , and 



(i) 



5. A 



Km ffl if^ooT^i = 7g SjA - Hence, sup K 7^ = 7^. Taking this into the account in ([U7j), 

we derive 



7^, A > sup 

*S,A.S» SA ^0 



— E (^ 1 ,...,s Jv) (Ai,...,A iV ] 



si,...,s. 



£s,A 



(C9) 



This proves equality (1Q2|) if set A is infinite and countable. 
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(b) Continuous case. Let, for a correlation scenario £s,a, each outcome set A n be infinite 
and uncountable. For positive integers K n > 1, introduce partitions 

Ufe^f^An, fc n = l,..,2*"} 



(CIO) 



of each set A n , such that 



7-)(2fcn 1) I | 7-)(2fcn) p<(kn) I 



fc n — 1 , . . . , 2 



VAT n e N. 



(Cll) 



For some partitions D^ i; D^ JVi of sets Ai, Ajy, respectively, consider collections ^cf^ ' ^ 
of real- valued measurable functions 

$£S*)(Ai.--a*):= X) ^:;:3^)( A i)-----x D c^(AAr), (ci2) 



specified by tuples A" := (K 1 , ...,K N ) and collections f3 K := {/3( s * 



(fei,...,fcjv) 

v) 



S n 1 , . . . , S^, 



fc n = 1, 2 Kn } of real numbers. Similarly to the derivation of (107|) . we have: 



> SU P 



^- E (V'(. 1 ,...^)(Ai,...,A^)) g ^ 



S1,...,S. 



> sup 



-j— E (^2 w) (Ai,...,A W : 



B~(J3 K , K) 



^S.A 



(K) 



where 



(C13) 



7 ^:=inf{e + + ne ± >o, e + -r = i, ^jy,^ 



Z/m; o Ihv 
S,A 



(si,...,s„) \ *i 



t~ p^S,A U ) / n (A:i) 
S ^(3i,...,3„) ^ATi 



x 



x ••• x D 



D™ x ... x 



K N 



(C14) 



* (si,...,s n ) y 



for each A" = (K u ...Kn), K n e N. From flC3l) . (041) and the special construction (Olj) 



of partitions Dk, it follows that jgp A < ls^ A < 7£ S a' ^ ^ — Li,...,Kn < Ln, and 

Substituting this into (Q3j) . we 



lim^,...,^^ 7^ 



(L) 

S,A 



7f s , A - Therefore, sup K 7^ s ; = 7^. 
prove equality (IC2p in case of uncountable sets A n . 

Coming back to notation £s,h ~ >* £p,M SA an d taking supremum of the left-hand and the 
right-hand sides of ()C2j) over all collections Ms a of POV measures, we prove relation f )52|) . 
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